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Abstract — This paper applies results on the robust stability 
of nonlinear quantum systems to a system consisting an optical 
cavity containing a saturated Kerr medium. The system is 
characterized by a Hamiltonian operator which contains a non- 
quadratic term involving a quartic function of the annihilation 
and creation operators. A saturated version of the Kerr non- 
linearity leads to a sector bounded nonlinearity which enables 
a quantum small gain theorem to be applied to this system in 
order to analyze its stability. Also, a non-quadratic version of 
a quantum Popov stability criterion is presented and applied 
to analyze the stability of this system. 

I. Introduction 

The use of Kerr media is commonly found in applications 
of nonlinear optics; e.g., see [1], [2]. A Kerr medium is 
characterized by a refractive index which increases with the 
intensity of light applied to the medium; e.g., see Section 
9.1.1. of [3]. Within the area of quantum optics, a Ken- 
medium is often characterized by a Hamiltonian operator 
which is a quartic function of the annihilation and creation 
operators; e.g., see Section 5.4 of [4]. This leads to a 
nonlinear quantum stochastic differential equation which 
contains a cubic nonlinearity [4]. In this paper, we apply 
some recent and new quantum robust stability analysis tools 
to analyze the stability of an optical cavity containing a Kerr 
medium. Such as system has been proposed as a method of 
generating squeezed light; see Chapter 9 of [3]. Squeezed 
light is an intrinsically quantum phenomenon which has 
potential applications in areas such as gravity wave detection, 
precision metrology and quantum computing [3], [5]. Also 
note that the quantum dynamics obtained in the case of a 
microwave resonator containing a Josephson junction can be 
used to approximate the case of a Kerr medium in a cavity; 
e.g., see [6]. Such a system is related to the system analyzed 
in [7]. 

The first method we will apply to analyze the robust 
stability of the system under consideration is the quantum 
small gain result presented in [8]. This result gives a suffi- 
cient condition for the robust stability of uncertain nonlinear 
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quantum systems in which the uncertainty is introduced 
by considering a non-quadratic perturbation to the system 
Hamiltonian operator Such a non-quadratic perturbation 
leads to a nonlinear quantum stochastic differential equation 
describing the system; e.g., [9]. This nonlinearity is required 
to satisfy a certain sector bound condition. Related results to 
the results of [8] can be found in [10], [11] and consider 
different classes of perturbations. Furthermore, the paper 
[12] introduces a quantum version of the Popov stability 
criterion (e.g., see [13] for the classical Popov stability 
criterion), which allows for quadratic perturbations to the 
system Hamiltonian. That is, the paper [12] considers un- 
certain quantum linear systems. In this paper we introduce 
a new version of the quantum Popov stability criterion 
which allows for non-quadratic perturbations to the system 
Hamiltonian and thus nonlinear uncertain quantum systems. 
As in [8], the nonlinearity is required to satisfy a certain 
sector bound condition. This result is applied to analyze 
the system consisting of an optical cavity containing a Ken- 
medium. 

For the quantum robust stability result introduced in [8] 
and in the new quantum Popov stability result introduced in 
this paper, the nominal quantum system is assumed to be a 
quantum linear system; e.g., see [14]-[18]. In addition, the 
nonlinearity is required to satisfy certain sector bound and 
smoothness conditions. However, for the standard quartic 
Hamiltonian model of a Ken medium, the resulting cubic 
nonlinearity will not satisfy the sector bound conditions for 
any finite sector. We overcome this difficulty by noting that 
any practical implementation of a Kerr medium will not 
be precisely modelled by a quartic Hamiltonian but rather 
will suffer from some saturation effects; e.g., see [19]. This 
allows us to model the Ken medium with a non-quadratic 
Hamiltonian such that the sector bound and smoothness 
conditions required in our quantum robust stability analysis 
results are satisfied. 

The remainder of the paper proceeds as follows. In Section 
Hn we define the general class of nonlinear uncertain non- 
linear quantum systems under consideration. In this section, 
we also also recall the main result of [8] and present a new 
Popov type stability result for this class of nonlinear quantum 
systems. In Section |III1 we analyze the system consisting 
of an optical cavity containing a saturated Kerr nonlinear- 
ity using the two quantum robust stability analysis results 
presented. In Section |IV] we present some conclusions. The 
proofs of all of the main results are given in the Appendix. 



II. Robust Stability of Uncertain Nonlinear 
Quantum Systems 

In this section, we describe the general class of quantum 
systems under consideration. As in the papers [8]-[10], [12], 
[20], we consider uncertain nonlinear open quantum systems 
defined by parameters {S,L,H) where S is the scattering 
matrix which is typically chosen as the identity matrix, L 
is the coupling operator and H is the system Hamiltonian 
operator which is assumed to be of the form 
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Here a is a vector of annihilation operators on the underlying 
Hilbert space and a# is the corresponding vector of creation 
operators. Also, AI E c2nx2n j^ ^ Hermitian matrix of the 
form 

A/i A/2 
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M* 



M* 



(2) 



and Ml = Afj* , A/2 = A/J. In the case vectors of operators, 
the notation ^ refers to the transpose of the vector of adjoint 
operators and in the case of matrices, this notation refers 
to the complex conjugate transpose of a matrix. In the case 
vectors of operators, the notation "^ refers to the vector of 
adjoint operators and in the case of complex matrices, this 
notation refers to the complex conjugate matrix. Also, the 
notation * denotes the adjoint of an operator The matrix A/ 
is assumed to be known and defines the nominal quadratic 
part of the system Hamiltonian. Furthermore, we assume 
the uncertain non-quadratic part of the system Hamiltonian 
/(z, z*) is defined by a formal power series of the form 



00 00 



/(z,z*) = ^5]5,,z^(z*)^ 
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which is assumed to converge in some suitable sense. Here 
Ski — S^f., and z is a known scalar operator defined by 



z = Eia + E2a* 
= [El E2] 
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The term /(z, z*) is referred to as the perturbation Hamil- 
tonian. It is assumed to be unknown but is contained within 
a known set which will be defined below. Two different sets 
of perturbations will be considered depending on the robust 
stability condition which is to be applied. 

We assume the coupling operator L is known and is of 
the form 



[ iVi N2 



where Ni e C"^" and N2 <E 



a 

n* 



(5) 



Also, we write 
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The annihilation and creation operators are assumed to 
satisfy the canonical commutation relations: 
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where J ~ 



= J 
e.g., see [18], [21], [22]. 
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To define the set of allowable perturbation Hamiltonians 
/(•), we first define the following formal partial derivatives: 
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Then for given constants 7 > 0, /3 > 0, (5i > 0, (^2 > 0, 
63 > 0, we consider the sector bound conditions 
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and the smoothness conditions 
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dzdz* dzdz 

Also, we consider the following upper and lower bounds on 
the perturbation Hamiltonian 



< f{z,z*) <Pzz* 



(14) 



Then we define two possible sets of perturbation Hamil- 
tonians >Vi and 'W2 as follows: 



/(•) of the form (|3]l such that 
conditions ( fTOl i and ( fT2] i are satisfied 
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/(•) of the form Q such that 
W2= { conditions O, (Hll, O and ([III ) . (16) 
are satisfied 

As in [8], [10], [12], we consider a notion of robust mean 
square stability. 

Definition 1: An uncertain open quantum system defined 
by iS,L,H) where H of the form Q, /(•) G W, and L 
of the form ^ is said to be robustly mean square stable if 
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Here 



ait) 
a*it) 



denotes the Heisenberg evolution of the 



vector of operators ^ ; e.g., see [20]. 

The following small gain condition is sufficient for the 
robust mean square stability of the nonlinear quantum system 
under consideration when /(•) G Wi: 

1) The matrix 

F = -iJM - -JNUN is Hurwitz; (18) 
2 

2) The transfer function 

Gis) = 2iE*T, isl - F)~^ J:JE^ 
satisfies the H°° norm bound 
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This result is given in the following theorem which is 
presented in [8]. 

Theorem 1: Consider an uncertain open nonlinear quan- 
tum system defined by (S*, L, H) such that H is of the form 
([ill, L is of the form (|5]l and /(•) G Wi. Furthermore, assume 
that the strict bounded real condition (fTSl l. (l20t is satisfied. 
Then the uncertain quantum system is robustly mean square 
stable. 

In the next section, we will apply this theorem to analyze 
the robust stability of a nonlinear quantum system corre- 
sponding to an optical cavity containing a Kerr medium. 
We also consider a new sufficient condition for robust mean 
square stability when /(•) G "¥^2, which is a nonlinear 
quantum version of the Popov stability criterion. This new 
condition is the existence of a constant 6 >Q, such that the 
matrix F defined in ( fTSl ) is Hurwitz and the transfer function 
Gis) defined in ( fT9l ) satisfies the strict positive real condition 

7 + (1 + eiuj)Giiuj) + (1 - eiuj)GiiLo)* > (21) 

for all ijj E [—00,00]. This result is given in the following 
theorem. 

Theorem 2: Consider an uncertain open nonlinear quan- 
tum system defined by (S*, L, H) such that H is of the form 
([T]i, L is of the form (|5]l and /(•) G >V2- Furthermore, assume 
that there exists a constant 6 > Q such that the matrix 
F defined in ( fTSl ) is Hurwitz and the frequency domain 
condition (ISTT i is satisfied. Then the uncertain quantum 
system is robustly mean square stable. 



In order to prove this theorem, we require the following 
definitions and lemmas. 

Lemma 1 (See Lemma 3.4 of [20].): Consider an open 
quantum system defined by iS, L,H) and suppose there ex- 
ists a non-negative self-adjoint operator V on the underlying 
Hilbert space such that 

- t[V,H] + h^V,L] + i[Lt,y]i + cV<X (22) 

where c > and A are real numbers. Then for any plant 
state, we have 

(Vit)) < e,-"* (V) + -, Vt > 0. 
c 
In the above lemma, [•, J denotes the commutator between 

two operators. In the case of a commutator between a scalar 

operator and a vector of operators, this notation denotes the 

corresponding vector of commutator operators. Also, Vit) 

denotes the Heisenberg evolution of the operator V and (•) 

denotes quantum expectation; e.g., see [20]. 

We will consider "Lyapunov" operators V of the form 
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where P G C^"^^" is a positive-definite Hermitian matrix 
of the form 

'Pi P2 ' 
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and 6 > Q. Hence, we consider a set of non-negative self- 
adjoint operators V defined as 

V of the form (|23ll such that P > is a 
Hermitian matrix of the form (l24l 

(25) 

Lemma 2: Given any positive definite matrix P of the 
form (EHi, then 
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which is a constant. 

Proof: The proof of this result follows via a straightforward 
but tedious calculation using @. D 

Lemma 3: With the variable z defined as in (|4|i and L 
defined as in (|5]), then 
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which is a constant vector. 
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In addition 
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which are constants. 

Proof: The proofs of these equations follows via straightfor- 
ward but tedious calculations using (|6]l. □ 
Lemma 4: Given any Hermitian matrix P of the form 
(l24l i. then the Hermitian operator 



Therefore using (|29] | and (l30l l. it follows that 
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which holds for any k > 0. Similarly 
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Now given any A: > and ^ > 0, let Hki = z^{z*Y and 
we have 
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and the constant /i is defined as in 
Proof: First, given any A: > 1 note that 
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Therefore, given any /(z, z*) E W2, 
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Hence using ( |28] |, the condition ( l27b is satisfied. D 

(30) Lemma 5: Given any /(z, z*) e W2 and L defined as in 
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a constant. Then, 
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Proof: In a similar fashion to the proof of Lemma |4] we 
write 
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It follows using ( l34l i that the condition ( |33] ) is satisfied. □ 

Lemma 6: Given a positive definite matrix P of the form 

(|24] |. a Hermitian matrix M of the form (|2]i, and L defined 
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Proof: The proof of these identities follows via straightfor- 
ward but tedious calculations using (|6]i. □ 
Lemma 7: Suppose z is defined as in (HJi and L is defined 
as in (|5]l. Then for any positive definite matrix P of the form 
(I24I 1 and any Hermitian matrix AI of the form (|2]i. 
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Proof: The proof of these equations follows via straightfor- 
ward but tedious calculations using (|6]l. □ 
Lemma 8: Given a complex row vector T = [Ti T2] . Then 
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Proof: The proof of this result follows via straightforward 
calculations. □ 

Proof of Theorem |2] If the conditions of the theorem are 
satisfied, then the transfer function ^ — (l+0s)G(s) is strictly 
positive real. However, this transfer function has a state space 
realization 
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E'^Tj. It now follows using the strict positive real lemma 
that the linear matrix inequality 
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will have a solution P > of the form (l24l i. This matrix 
P defines a corresponding Lyapunov operator V^ e P as in 
(I23I 1. Furthermore, it is straightforward to verify that HG + 
G^H^ = 0. Hence, using Schur complements, it follows 
from (|44]i that 

PF + F^P 

^ {PG + H^ + OF^H^) [G^'P + H + OHF) < 0. 
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Now using Lemma [T] and Lemma [8] we have 
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where 
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F = ~iJM - \JNUN, G = 2iJj:E^ and H = P#E. 
Also, it follows from (|45j that M < 0. 

We now write V E V as 
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where 



V 



a' a 



]P 



a 

n* 



Also, we define 



H 



a' a 



M 



a 



Hence (l46b can be re-written as 



+ ^[ +l{Lnz,L] + [L^,z]L)+z 

i[z,V-eH] 
ff (Lt[z,i] + [Lt^^]L)+; 



a 



M 



Tr ( P JA^t 



a 

n* 



I 





NJ 



(47) 



Now, it follows from Lemma |4] and Lemma |5] that since V — dH is self-adjoint. Therefore, 

-. -. / i[z.V-0H] + ^(L^lz,L] + [L\z]L) \ 

_2 2_ y +^{wi-z*/-/) J 

= -i[V + efiz,z*),H + f{z,z*)] , ^[zy-eH]+i{L^[zM+[L'.AL) \* 

+\L^[v + efiz^z*).L] + \[L\v + efiz.z*)]L -[- +^i^,'lz*MT ) 

= -z[V,H]-ze[fiz,z*),H]~^[VJ{z,z*)] 1/.^ - - e ,.^ ^^ ,. ^ . 

1 1 = -[i[z,V-0H] + -{U[z,L] + [Llz]L, 

~ie[f{z,z*)J{z,z*)] + -L^[V,L] + -[L\V]L l\ 2 

4LHfiz,znM^?,iLKfiz,zn]L x(.[.^^-.i^] + ^(Lt[.,,H[^^.1x] 



-i[V, H] - i[V - 9H, f{z, z*)] _ f _ \[^' "'^ - ^^j 



^ ^ +^{wi~z/^) (wi~z-/^) 

-z[V,H] + -L^[V,L] + -[L\V]L 

-i[V - 9H, z]wi + iwl[z*,V - 9H] and hence 

i r r- „ - -,n i 



-- [z, [V - dH, z]] W2 + -wl [z, [V - OH, z]] 



+ -wl {L^[z*,L] + [L\z*]L) +wl (i[z*,V ^eH] + ^ {l''[z*, L] + [L\ z*]L) 

--[[L\z]L,z]w2 < -(■i[z,V-0H] + -(L''\z,L] + \L\z]L] 

Q 7 V 2 ^ ^ 

~-w*,[[L\z]L,z]* ,,Uz*,V-6H] + Ul\z\L] + [L\z* 

9 V 2 

+ - [z*, [L\z]L\ W3 +^ (y,^ _ ^*/7)* K - 277) 



+ -ti;* [z*, [L\z]L] U[z,V-eH] + - {L^[z,L] + [L\z]L) ) z* h 



1 .4,.-. ., 1 

-I 



{V,H\ + -it[^,L] + -[Lt,y]L ^^ Uz*,y - 0^] + ^ {L\z\L\ + [Lt,z*]L) ) /7 



2"'H +e[[in2]L,2]* 

-r\B{z\\L\z\L\w3 



7 



Furthermore, 



2 

--0w;3*[z^[L^z]L]^ (48) o < ({B\^L\z\L,z\^i{z,^-dE,z^)* ^W2 

x((4[i^z]L,z]+z[z,[\/-0iJ,z]])*+u.2) 
= (0[[Lt,z]L,z]+t[z,[i^-0iJ,z]]) 
X (0[[Lt,z]L,z]+z[z,[y-0i?,z]])* 
+i;;2*(0[[Lt,z]L,z]+z[z,[y-0i/,z]])* 
"^ -BB,z\* =z*{V-9H)-(y-0H)z* = [z*,V-9H] + {9 [[L\z]L, z] + i[z,[V - 9H, z]]) W2 + W2W2 



Also, 



and hence 



(61 [[l\ z\l, z] + i [z, [V - en, z]] ) W2 

wl{e[[L\z]L,z\* -i[z,[V-6H,z]X 
< {d[[L\z]L,z]+i[z,[V~eH,z]]) 

X {0 [[L\z]L, z] + 1 [z, [V - eii, z]])* 

+W2W2. (50) 



Similarly 



< i^e[z*,[L\z]L]* ~W3^ 

x(^9[z*,[L\z]L]*~W3^ 
= 6|2 [z^[L^z]L] [z*,[Lt,z]L]* 

-w;e [z*, [L\z]L] * - e [z*, [l\z]l] W3 + u.>3 



and hence 



e [z*, [l\z]l] w3 + w^e [z*, [l\ z]l]* 

< 9^z*,[L^,z]L][z*,[L\z]L]* +w;w3. (51) 
Substituting (|49]), (|50| and (|5T]l into (gS]), it follows that 

-^[V,H] + h^[V,L] + ^[L\V]L 
< -i[V,H] + ^L^[V,L] + ^[L\V]L 

+ - (i[z, V-eH] + ^ {L^z, L] + [L\ z]L) + z\ 

X (t[z*, V-eH] + ^ {L^[z*,L] + [L\z*]L) + z* 

zz 
+7 {wi - z*/7)* {wi - z*/7) 

7 

+ ]-{e[[L\z]L,z]+i[z,[V-eH,z]]) 



X 



e[[L\z]L,z]+f[zAV~OH,z]] 



using (HB, (O, O, and (O. Then it follows from (|47] | 
that 



-i[V,H] + \l^[V,L] + \[L\V]L 



< 



a 



M 



a 

n* 



Tr I P.JN'^ 

fl2 



/ 




NJ 

[z\[L\z]L\[z\[L\z]L\ 



{6[[L\z\L,z\ +i[z,[V~eH,z\]) 
[[L\z]L,z\+i[z,[V-eH,z]]Y 






Since M > 0, it follows using (fT4l) that there exists a 
constant c > such that 



-t[V, H] + Il^ [V, L] + ^[L\V]L + cV 

< -i[V, H] + iit [V, L] + ^[L\V]+c{V + e(3z*z) 



< Tr PJN"^ 



I 




NJ 



1 * 

+ -^2^2 

+ ^-^[z*,[L\z]L][z*,[L\z]lY + ^w;ws 
< -t[V,H] + h^[V,L] + ^[L\V]L 

+ - (i[z,V - eS] + - {L''[z, L] + [L^ , z]L) +z 

X I t[z*, V-eH] + ^ {L^[z*,L] + [Llz*]L) + z* 



+ -[z\[L\z]L\[z\[L\z]lX 

+ ]^{6[[L\z]L,z\+i[z,[V~eH,z]]) 
X {6[[L\z]L,z]+i[z,[V-eH,z]]y 
+^17+1 + 1 (52) 



+ "-[z*,[L\z]L][z*,[L^,z]L]* 

+ ^{d[[L\z]L,z]+i[z,[V~9H,z]]) 

X {e [[L^,z]L, z] + i [z, [V - en, z]])* 

+'5i7+y + y. 



That is. 



i[V, H] + ^L^V, L] + i [L\ V]L + cV <\ 



where 

A 



= Tr PJN^ 



I 




NJ 



-Ue[[L^,z]L,z]+i[z,[V 



]L][z*,[L\z]L] 

-9H,z]]) 
dH,z]]y 



-Si'y 



[[L^,z]L,z]+i[z,[V 
2 2 



= Tr 



(pjN^ 



32 



/ 




NJ 



— E* jj:n^ n*j:je^ ejj:n^ NY.JE^ 

2 



-OEJN^NT.JE'^ 
+iEY.J (P - f A/) JE'^ 

-OEJN'^NY.JE'^ 
-iEJ:j {P - f A/) JE^ 

2 2 



-(5i7 + TT + TT 



> 



using Lemma |2] and Lemma [3] Therefore, it follows from 
Lemma [U and P > that 



a{t) 
a* it) 



a{t) 
a* it) 



< 



a(0) 
(0) 



it 



i# 



«(0) 
a#(0) 



\„^aAP + l3E^E] 



^min [P\ 



A 



C'^min [-* J 



vt>o. 



(53) 



Hence, the condition 



( fTTb is satisfied with ci 



> 0, C2 = c > and C3 = 



A 

cA™.„[P] 



> 0. 



Observation 1: Note that the SPR condition ( |2TI ) can be 
re-written as 

^+7^e[G(ia;)]-6lwZm[G(iw)] > (54) 

for all id e [—03,00]. The condition ( l54l i. can be tested 
graphically by producing a plot of ojIm[G{iuj)] versus 
TZe[G{iuj)] with ut € [— cxj,cxj] as a parameter Such a 
parametric plot is referred to as the Popov plot; e.g., see 
[13]. Then, the condition ( l54l i. will be satisfied if and only 
if the Popov plot lies below the straight line of slope | and 
with a--axis intercepts 



-^; see Figure [T] 



III. Analysis of an optical cavity containing a 
Kerr medium 

The system under consideration consists of an optical 
cavity containing a Kerr medium. The optical cavity is 
made from two mirrors, one of which is partially reflecting 
and one of which is fully reflecting. The cavity is driven 
by a laser beam directed at the partially reflecting mirror 
The corresponding reflected beam is then measured using 



allowable region 



ujTm[G{iLj)] 



slope = ^ 




ne[G{iuj)] 



Fig. 1. Allowable region for the Popov plot. 



a detector The Kerr medium within the cavity can be 
constructed from a suitable nonlinear optical crystal; e.g., 
see [3]. This system is illustrated in Figure |2] 



Partially reflecting 
miiror 




Fully reflecting 
minor 



Kerr medium 



Fig. 2. Schematic diagram of an optical cavity containing a Ken' medium. 

A Standard {S, L, H) model for an optical cavity contain- 
ing a Kerr medium is as follows: 



S ^ I, H = (a*) a^, L = y/Jia; 



(55) 



e.g., see [4]. We first attempt to apply the results of Theorem 
[T] and Theorem |2] to this quantum system. Hence, we let 



and 



A/ = 



f{z,z*)^z^iz*f 



where z = a*. This defines a nonlinear quantum system 
of the form considered in Theorem [1] and Theorem |2] with 

A/i = 0, A/2 = 0, iVi = y/J^, N2 ^ 0, El ^ 0, E2 = 1. 



We now investigate whether this function /(•) satisfies the 
conditions ([Hi, ^, O, ^, and ([T3]l- Now, 



dz 

d^f{z,z*) 

dz^ 

d^f{z,z*) 

dzdz* 



2z{z*y 

2{z*f- 



Azz* 



Also, the sector condition (fTTT i can be rewritten as 

df{z,z*rdf{z,z*) 



< 



dz dz 

df{z,z*) 



dz 



z + z 



df[z,z*) 



dz 



7(5i 



which is satisfied for 7 = since 



df{z,z*Y 



z + z 



df{z,z*) 



Az'^{z*Y>Q. 



dz dz 

However this condition is not satisfied for any finite value of 
7 > 0. Also, conditions ([TOl), (O, ([Tl]i, ^^, ([T3) are not 
satisfied. 

In order to overcome this difficulty, we note that any 
physical realization of a Kerr nonlinearity will not be exactly 
described by the model (ISST i but rather will exhibit some 
saturation of the Kerr effect; e.g., see [19]. In order to 
represent this effect, we will assume that the true function 
/(•) describing the Hamiltonian of the Kerr medium is such 
that its Taylor series expansion Q satisfies 5o.fe = Si^k = 
for all k = 0, 1, . . . and ^2.2 = 1- That is, the first non- 
zero term in the Taylor series expansion corresponds to the 
standard Kerr Hamiltonian given in dSST l. Furthermore, we 
assume that the function /(•) is such that the conditions ( fTOl i. 
(fTTT i. (fT2l i. ( fT3l l. ( fT4l i are all satisfied for suitable values of 
the constants j > 0, P > 0, Si > 0, 62 > 0, S3 > 0. Here 
the quantity - will be proportional to the saturation limit for 
the Kerr nonlinearity. Thus, under these assumptions, we can 
assume /(•) e Wi and /(•) e ^2- 

^ 



This system has F = 





2j 



which is Hurwitz 



for all K > and G(s) = --rV. A magnitude Bode plot of 
this transfer function, is shown in Figure |3] for the case of 
K = 2. In this case, we obtain ||G'(s)||oo = 2 and in general 

l|G(5)|U = -. 

K 

Thus, applying Theorem [T] to this system, we can guarantee 
that the system is mean square stable provided 

4 



K > 



7 



(56) 



We now apply our new result Theorem |2] to further 
analyze the stability of the system. We first choose k = 2 
and construct the Popov plot corresponding to the transfer 
function G{s) as discussed in Observation [T] For a value 
of 6* = 1, this plot, along with the corresponding allowable 
region corresponding to 7 = 0.1, is shown in Figure|4] From 
this figure it can be seen that the Popov plot lies in the 



Frequency (rad/sec) 



Fig. 3. Magnitude Bode plot of G{s) for the case k = 2. 



allowable region and hence, it follows from Theorem |2] and 
Observation [T] that this system will be mean square stable 
for K = 2 and 7 = 0.1. In fact, it follows from this plot 
that the frequency domain condition ( 1211 1 will be satisfied 
for all 7 > 0. This condition is clearly less restrictive 
than the condition ( l56b obtained by applying Theorem [T] 
Furthermore, we can construct the Popov plot of the system 
for different values of k > as shown in Figure [5] From 
these plots, we can see that for a suitable value of 6 — - > 0, 
the frequency domain condition (ISTT i will be satisfied for 
all 7 > and all k > 0. Thus, using Theorem |2] we can 
conclude that the optical cavity containing a saturated Kerr 
medium is in fact mean square stable for all 7 > and 

K > 0. 



,3 1 



E 0.5 

H 




■R.a[G{iuj)] 



Fig. 4. Popov plot for the Kerr nonlinearity system with k 
7 = 0.1. 



2 and 



IV. Conclusions 

In this paper, we have introduced a new nonlinear quantum 
Popov stability criterion and applied it to the robust stability 




Fig. 5. Popov plot for the Ken' nonlineaiity system with different values 
of /i > 0. 



analysis of a nonlinear quantum system consisting of an op- 
tical cavity containing a Kerr medium. We have also applied 
an existing quantum small gain theorem to the analysis of this 
system. By choosing a model which represents a saturating 
Kerr medium, both approaches to robust stability analysis 
were applicable to this system. Furthermore both approaches 
were able to verify the robust mean square stability of this 
system for some range of parameter values. However, the 
quantum small gain theorem approach was found to be more 
conservative than the quantum Popov criterion approach in 
that it could only verify robust mean square stability for 
a restricted range of parameters. In contrast, the quantum 
Popov approach was able to verify the robust mean square 
stability of the system for all positive values of the system 
parameters. 
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